syms x y
p = expand((x+1)^3)
q = expand(sin(x+y))
